A Simulation Study of the Lithium Ion Transport Mechanism in Ternary Polymer 
Electrolytes - The Critical Role of the Segmental Mobility 
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The lithium transport mechanism in ternary polymer electrolytes, consisting of poly(ethylene 
oxide) (PEO), lithium-bis(trifluoromethane)sulfonimide (LiTFSI) and the ionic liquid N-methyl-N- 
propylpyrrolidinium bis(trifluoromethane)sulfonimide (PYR13TFSI) is analyzed vi molecular dy- 
namics (MD) simulations. Starting from the classical, binary electrolyte PE02oLiTFSI, we focus 
on two different strategies by which the ternary electrolytes can be devised, namely by (a) adding 
the ionic liquid to the binary system, and (b) substituting the PEO chains in PE02oLiTFSl by 
the ionic liquid. In order to elucidate in how far the microscopic lithium transport differs between 
these two electrolyte classes, we employ an analytical. Rouse-based cation transport model^, which 
has originally been devised for the binary systems. Within this framework, the cation transport 
is characterized via three different mechanisms, each quantified by an individual time scale. Our 
analysis demonstrates that this model is also applicable to the ternary electrolytes, allowing us to 
express the effect of the ionic liquid in terms of the time scales of the transport model. Whereas the 
addition of the ionic liquid to PE02oLiTFSI plasticizes the polymer network and thus also increases 
the lithium ion mobility, no significant effect can be observed when substituting the PEO chains 
by the ionic liquid. This is due to the fact that in the latter case the amount of free, mobile ether 
oxygens reduces, since more lithium ions coordinate to the PEO backbone, which compensates with 
the plasticizing effect. Thus, the segmental mobility displays a decisive role in polymer electrolytes. 
In total, our findings are agreement with recent experimental observational. 

PACS numbers: 



I. INTRODUCTION 



Polymer electrolytes - typically consisting of an amor- 
phous polymer matrix and a salt dissolved in it - 
are promising candidates for many technological devices 
such as batteries or fuel cells. They were first discov- 
ered by Wrigth in 19733-, who investigated salt-doped 
poly(ethylene oxide) (PEO) melts. Some years later their 
potential for electrochemical applications was discovered 
by Armand^, since they combine many properties of con- 
ventional liquid and ceramic electrolytes due to their in- 
herently stable but flexible nature^i^. However, at ambi- 
ent temperatures the conductivity of many of these ma- 
terials is still too low for an efficient technological use. 

Several attempts have been made in the past to over- 
come this deficiency. Common modifications of the clas- 
sical polymer-salt systems are the addition of nanoparti- 
ciesX or organic solvent molecules^"—. However, in case of 
the latter, the resulting electrolytes suffer from the high 
volatility and thus fiammability of the solvent as well 
as its reaction with lithium metal electrodes*!^. Here, 
Passerini et aZ.— liir— demonstrated that the use of an 
ionic liquid instead of a conventional solvent has sev- 
eral advantages: The ionic liquids are non-volatile, non- 
flammablei^ and exhibit a wide electrochemical stability 
windowi^. In this way, improved polymer electrolytes 
can be designed, which show an increased conductivity 
combined with inherent stability^i^ and are thus ideal to 
create light-weighted but powerful batterie o^^'^'^ . 



However, it is not yet fully understood in how far the 
lithium transport mechanism in these materials changes 
relative to the conventional polymer electrolytes. For in- 
stance, it was speculatedii that the lithium ions become 
progressively coordinated by the anions from the ionic 
liquid and are thus decoupled from the rather slow PEO 
chains. Alternatively, one might also expect that the IL 
enhances the PEO dynamics and serves as a plasticizer 
in this way, which is a common observation when adding 
low-molecular solvents to PEO-salt systems^""—. In this 
contribution, we utilize molecular dynamics (MD) simu- 
lations to address these issues, and interpret the changes 
in the microscopic lithium transport mechanism within 
the framework of an analytical model^'--^, which has orig- 
inally been devised for binary polymer electrolytes. 

With few exceptionsi^r— , the archetypal PEO-based 
electrolytes have been most intensely studied via MD 
simulationsi^i^ir— . Here, essentially three different 
cation transport mechanisms have been observed: dif- 
fusion of the cation along the polymer chain, cooperative 
motion of the cation with the polymer chain and cationic 
transfer between different polymer chains. While all of 
these mechanisms have already been observed in early 
studies^i2ij only relatively recent work has resulted in a 
quantitative characterization of the individual contribu- 
tions to the overall cationic difFusioo^ '^^i^^ . This is due to 
the fact that several ten to hundred nanoseconds are re- 
quired in order to fully quantify the relevant microscopic 
transport processes. Moreover, an appropriate force field 
is important in order to gain reliable results. Here, espe- 
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Figure 1: Sketch depicting the three different cation trans- 
port mechanisms in PEO-salt electrolytes. Each mechanism 
is characterized by a specific time scale. 



daily the incorporation of polarization promises a signif- 
icant improvement in the description of the ion dynam- 

Based on insights from simulations of a PEO/LiTFSI 
electrolytes^, Borodin et al. have proposed a transport 
model based on the three mechanisms that was consistent 
with the numerical MD data. Moreover, the model was 
able to predict the experimentally observed lithium diffu- 
sion coefficient. To this purpose, the diffusion coefficient 
for the motion of the ion along the polymer chain, the 
intersegmental jump rate (i. e. the rate of hopping pro- 
cesses between segments of two different PEO chains or 
between two remote segments of one PEO chain) as well 
as the characteristic form of the segmental MSD between 
two jumps, i.e. MSD(i) = a{c,T)t"-^ (with a being a 
function of the salt concentration c and the temperature 
T), were extracted from the simulations. These three in- 
gredients were combined in a microscopic model similar 
in spirit to the Dynamic Bond Percolation (DBP) model 
of Ratner— . By utilizing a Monte- Carlo scheme for the 
motion along the chain and the intersegmental hopping 
process, the lithium diffusion coefficient was obtained nu- 
merically from this model. 

Motivated by the success of this model, Maitra and 
Heuer developed an analytical framework to describe the 
cation transport in a PEO/LiBF^ electrolyte^, which is 
based on both the Rouse modeP^ as well as the DBP 
model^^. Here, the purely empirical description of the 
segmental dynamics, MSD(t) — a{c,T) t'^-^ , as well as 
the numerical evaluation of the diffusion along the chain 
and the intersegmental hopping were replaced by ana- 
lytical expressions. Each of the three transport mecha- 
nisms was characterized by a specific time scale (see also 
sketch in Figure [1]): 1. Diffusion along the PEO chain, 
which can be interpreted as an effective one-dimensional 
random walk along the curvilinear path of the polymer 
chain. The time scale the ion needs to explore the en- 
tire PEO chain is denoted by ti. 2. Segmental motion of 
the PEO chain, which can be separated into the ccntcr- 
of-mass motion and the internal dynamics. The internal 
dynamics can be described by the Rouse model with an 
effective Rouse time T2 characterizing the motion of the 
bound PEO segments. 3. Intersegmental transfer of the 
cation from one PEO chain to another takes place. This 



motion can be quantified by the average residence time 
T3 at a given chain. It was showni that this mechanism 
can be viewed as a renewal process within the frame- 
work of the DBP model"^" , since the dynamics of a given 
lithium ion becomes independent of its past after being 
transferred to another PEO chain. Therefore, this con- 
tribution is also vital to the long-range cation transport 
and a significant conductivity on the macroscopic scale. 

Since the diffusion along the chain is also expressed 
within the framework of the Rouse model, correlations 
with the polymer dynamics are implicitly contained in 
this model, which have been neglected by Borodin et 
al.—. Another advantage is that the lithium diffusion 
coefficient can easily be calculated as a function of the 
polymer chain length N , which immediately follows from 
the iV-dependence of the three underlying time scales. 

Of course, as motivated above, the situation is ex- 
pected to change when other, more complex elec- 
trolytes than the binary PEO-salt mixtures are in- 
vestigated. In this work, we aim to generalize the 
transport model to ternary polymer electrolytes con- 
sisting of PEO/LiTFSI and the ionic liquid A^-methyl- 
A'^propylpyrrolidinium bis (trifiuoromethane) sulfonimide 
(PYR13TFSI, in the following also abbreviated as IL) 
in order to unravel the microscopic origin of the experi- 
mental observations by Passerini et al.^^. In particular, 
we focus on PE02oLiTFSI, to which a variable amount x 
of PYR13TFSI is added. However, it is a priori unclear if 
the lithium transport mechanism changes only quantita- 
tively - characterized by different values for ri , T2 and T3 
- or if it even changes on a qualitative level. Certainly, 
the transport mechanism will be altered at some point 
when PEO is successively substituted by PYR13TFSI 
molecules. In order to investigate such a crossover, addi- 
tional simulations with constant lithium volume concen- 
tration were performed, in which the PEO chains were 
successively substituted by PYR13TFSI molecules under 
the constraint of a constant volume of the simulation 
box. Figure [5] shows snapshots from the simulations for 
the two different scenarios. Since the transport mecha- 
nism in the binary PE02oLiTFSI electrolyte has already 
been thoroughly exploredi^^, this system serves as a 
reference substance. 

In what follows, the experimentally investigated elec- 
trolytesii with a constant ratio of EOs to lithium ions, 
namely EO : Li = 20 : 1, and a variable amount 
X of PYR13TFSI relative to LiTFSI will be denoted 
as PEOzoLiTFSI ■ x PYR13TFSI. The theoretically 
motivated systems to study the transition between bi- 
nary PEO/LiTFSI and PYRi3TFSI/LiTFSI electrolytes 
will be abbreviated as PEOso-a^LiTFSI- x PYR13TFSI. 
Here, a is the ratio of the respective partial molar vol- 
umes of PEO and PYR13TFSI and will be determined 
subsequently. 

For reasons of simplicity, we also abbreviate PEO as 'P' 
(i. e. polymer or PEO) and LiTFSI as 'S' (i. e. salt) in the 
following, leading to the short-hand notations P20S • x IL 
and P2o-cta;S • X IL for the two distinct classes of ternary 
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Figure 2: Snapshots from the simulation showing (a) the addition of IL to the binary polymer electrolyte PE02oLiTFSI and 
(b) the substitution of PEO chains by IL molecules from the same electrolyte. PEO chains are shown in red, lithium ions in 
green, whereas all other ions (PYR^g and TFSI~) are invisible. 



electrolytes. 



II. MD SIMULATION METHODOLOGY 

The MD simulations were performed with a modified 
version of the AMBER Molecular Dynamics Program 
Package^, allowing us to employ a specifically designed 
force field for PEO/LiTFSI^^-^s and PYRigTFSI^^. 

The initial configurations were created by randomly 
placing the individual molecules on a simple cubic lattice 
with a lattice constant of 36 A irrespective of their type, 
mixing the system in this way. In case of PEO, the chains 
already had coiled conformations as under melt condi- 
tions, thereby circumventing an expensive equilibration 
over several Rouse times. 

The binary PEO/LiTFSI electrolyte contained 10 PEO 
chains with TV = 54 monomers each as well as 27 
LiTFSI molecules. In case of the addition of IL, i.e. 
P20S • X IL, the simulation cell also contained 18 {x — 
0.66) or 87 {x = 3.24) PYR1 3TFSI molecules (the lat- 
ter is shown in Figure 2(a) I. In case of the substi- 
tution of PEO (P2o-Q2:S • X IL) and thus a constant 
lithium volume concentration, the partial molar volumes 
of PEO and PYR13TFSI must be known. Through sev- 
eral short equilibration runs of the binary PEO/LiTFSI 
and PYRiaTFSI/LiTFSI electrolytes, we found that 
both the molar volumes of PEO and PYR13TFSI are 
independent of the lithium mole fraction (at least for 
the investigated concentration range), and that approx- 
imately a = 7.14 PEO monomers occupy the same vol- 
ume as one PYR13TFSI ion pair. Based on these val- 
ues, we created the following systems: 8 PEO chains 
with 15 PYR13TFSI molecules (PigS • 0.556 IL), 6 PEO 



chains with 30 PYR13TFSI molecules (P12S • 1.111 IL) 
and 4 PEO chains wi th 45 PYR13TFSI molecules (PsS • 
1.667 IL, see Figure 2(b) I, all of them containing, as 



above, 27 LiTFSI ion pairs. Neat PEO (16 chains with 
N = 54) and the pure PYRigTFSI/LiTFSI electrolyte 
(76 PYR13TFSI as weU as 27 LiTFSI ion pairs) serve as 
reference within this context. 

The systems were equilibrated in the NpT ensemble for 
70 — 80 ns using the PME technique^. Afterwards, pro- 
duction runs with a total length of 200 ns were performed 
in the NVT ensemble. A time step of 1 fs was used in 
all simulations to propagate the systems. The tempera- 
ture was maintained by the Berendsen thermostat)^ to 
a reference temperature of To — 423 K. All bonds in- 
volving hydrogen were constrained by the SHAKE algo- 
rithm^^. The induceable point dipoles were integrated by 
the Car-Parrinello-like scheme^i, and the charge-dipole 
interactions between atoms separated by three bonds (1- 
4 interactions) were scaled by a factor of 0.2. For all 
other nonbonded contributions, the full 1-4 interaction 
was taken into account. Dipole-dipole interactions were 
damped using the Thole screening with ct — QA^^ . 

One indeed observes from Table U that the equilib- 
rium values for the box lengths are nearly the same 
for all P2o-Q2;S • X IL systems. Slight deviations for 
LiTFSI • 2.815 PYR13TFSI may resuh from rounded 
molecule numbers in the simulation. 



III. STRUCTURAL PROPERTIES 

In order to quantify the local structure around the 
lithium ions, radial distribution functions (denoted as 
g{r) in the following) have been computed for the atom 
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system 


Lbox [A] 1 PEO [%] 


2 PEO [%] 


6g [A^] 


[A^] 


(R^) [A^] 


PEO 




40.11 






10.39 


1904 ± 43 


306 ± 35 


P20S 




35.96 


47.2 


52.7 


9.71 


1573 ± 112 


260 ± 42 


P20S • X IL 


P20S 


• 0.66 IL 


38.27 


53.0 


47.0 


9.66 


1654 ± 121 


272 ± 38 


P20S 


■ 3.24 IL 


45.29 


75.8 


24.0 


9.58 


1498 ± 57 


249 ± 37 








P20- 


-axS ■ X IL 








PieS 


• 0.556 IL 


35.93 


55.4 


44.5 


9.54 


1767 ± 271 


267 ± 50 


P12S 


• 1.111 IL 


35.96 


64.2 


35.4 


9.34 


1349 ± 150 


227 ± 43 


PgS- 


1.667 IL 


36.00 


83.7 


13.8 


9.09 


1840 ± 430 


280 ± 80 


PoS- 


2.815 IL 


36.33 













Table I: Length of the simulation box and ratio of lithium ions coordinating to one or two PEO chains. Mean squared chemical 
bond length 69, mean squared end-to-end vector (Re) and mean squared radius of gyration (Rg) of the PEO chains. 



pairs Li — OpEO as well as Li — Otfsi (see supporting 
information). Both coordination types exhibit a sharp 
peak around 2 A corresponding to the first coordination 
shell, which is in good agreement with neutron diffraction 
experiments^ and quantum chemistry calculations^Siii. 
For Li"*" — OpEO, the first coordination sphere directly 
crosses over into a second small peak around 3.3 A. At 
larger distances no significant structural arrangement can 
be found. For the coordination of Li"*" — Otfsij peaks be- 
come noticeable also at larger distances around 6 A and 
8 A, thus demonstrating that long-ranged correlations 
are present as also observed for a PE02oLiI electrolyte^. 
The same observation can be made for the PEO-free elec- 
trolyte PqS ■ 2.815 IL, as also reported in previous MD 
studies^. The pair correlation function of the PYR13 
cations and TFSI anions (i. e. Npyri3 ~ Otfsi) exhibits 
only a weak first coordination peak between 3.6 and 
5.0 A. 

When successively adding the IL, one observes that 
the peak positions of both Li+ — Opeo and Li^ — Otfsi 
remain the same for all electrolytes. For P20S • x IL, 
the EO coordination numbers extracted from the inte- 
gral over the first shell increase slightly, partly as a re- 
sult of the lower fraction of lithium ions coordinating to 
two PEO chains (see discussion below), for which the 
EO coordination number is lower due to steric effects. 
In P2o-axS • X IL, the absolute EO coordination num- 
ber (as determined from the integral over g{r)) decreases 
with increasing IL content. This trend can be explained 
by the decreasing number of possible coordination sites, 
since the PEO concentration reduces. Naturally, the co- 
ordination number of the TFSI oxygens increases with 
increasing IL content in both types of electrolytes. 

Since the positions of the peaks do not change with 
composition, the same criterion to define temporary 
lithium bonds will be used for the subsequent analysis. 
We consider a EO and a Li"'" as bound if their distance is 
not larger than 3.0 A. In analogy, we consider a Li"*" and 
a TFSI oxygen as bound if their distance is not larger 
than 2.7 A. 

Figure |3] shows the probability distribution functions 



p{n) to find a lithium ion with n EOs or TFSI oxy- 
gens in its first coordination shell. One observes that 
the coordination numbers are very similar in all systems 
and only change slightly with the IL concentration. The 
predominant lithium coordination consists of 4 — 5 EOs, 
which is in good agreement with experimental dsXcS^ and 
quantum chemistry calculations^!^. In those complexes 
where the 4—5 EOs originate from a single PEO molecule, 
the polymer chain wraps helically around the cation. For 
complexes involving two PEO chains, typically 2 — 3 EOs 
from each chain coordinate to the ion. Coordinations by 
TFSI oxygens are rather rare, and in most cases the anion 
coordinates only briefly to the lithium ion (see also snap- 
shots in Figure O where all ions are in the vicinity of a 
PEO chain). The only exceptions are PgS • 1.667 IL and, 
though less pronounced, P12S • 1.111 IL, for which also 
an EO coordination number of n = 2 becomes noticeable. 
Simultaneously, the lithium ions more likely coordinate 
to 1 — 2 TFSI oxygens, while some are even coordinated 
to TFSI only (0.4 % for P12S • 1.111 IL and 2.5 % for 
PsS • 1.667 IL). As in pure IL/LiTFSI, the prevalent co- 
ordination number in the latter scenario is about 3 — 4, 
originating mostly from different anions. This specific co- 
ordination has also been observed previously in MD sim- 
ulations^'^, although experimental work emphasizes that 
also the Li(TFSI)2 complex is important^. 

Table H] summarizes the percentages of lithium ions co- 
ordinating to one or two PEO chains. Coordinations to 
three PEO chains were rarely observed and had a very 
brief life time of a few picoseconds only. Therefore, these 
events were neglected for the subsequent analysis. For 
the pure polymer electrolyte, the fractions of complexes 
involving one and two chains are nearly equal. Simi- 
lar binding energies for both coordination types have 
also been found in quantum chemistry calculations'*—. 
The fraction of lithium ions coordinating to one PEO 
molecule increases with the IL content for both types of 
electrolytes. This is a consequence of the reduced PEO 
concentration, as it becomes less likely that a lithium ion 
encounters a second chain in the semidilute case. 

Interestingly, in P20S ■ x IL, the fraction of both Li"*" 
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Figure 3: Probability p{n) to find a certain coordination number n of EOs (irrespective if tiie ion is tied to one or two PEO 
chains) or TFSI oxygens in (a) P20S • x IL and (b) P20-aiS • x IL. 



coordinating to TFSI and Li"*" coordinating to one PEO 
chain increases linearly with the IL amount x, thus indi- 
cating that the changes in the coordination sphere are 
purely statistical. For P2o-a2:S • x IL however, these 
trends are stronger than linear, which can be attributed 
to the fact that the PEO chains are successively removed 
from the system, and the remaining PEO chains become 
more crowded in this way. Thus, the different coordina- 
tion shell can not solely be explained as a simple dilution 
effect. 

Due to the helical coordination structure of the PEO 
backbone, the local polymer structure changes, and the 
conformational phase space of the chain is reduced. 
Moreover, in case of the ternary electrolytes, the ad- 
ditional IL molecules dilute the PEO molecules, thus 
inducing a crossover from a polymer melt to a semidi- 
lute solution, which may also alter the equilibrium con- 
formation of the polymer chains^^. Table |T] summa- 
rizes the mean squared distance 6q between two chemical 
monomers and the mean squared end-to-end vector (R^). 
Due to the crown-ether-like coordination of the PEO 
backbone, &q is smaller in all lithium-containing systems. 
For P2o-Q2:S • X IL, this trend becomes more pronounced 
with increasing IL concentration. A contraction of the 
polymer chain can also be observed for P20S • x IL from 
(R^), which also decreases, whereas for P2o-axS • x IL 
no clear predictions can be made within the error bars. 
When determining the ratio of (R^) and the radius of gy- 
ration (Rg), one finds values close to the ideal ratio of six 
for a Gaussian chain^^ for all electrolytes, again indicat- 
ing that Gaussian chain properties are fulfilled to a good 
approximation on a global scale. Also for the scaling of 
the Rouse mode amplitudes (X^), we observe only slight 
deviations from the respective curve for pure PEO (see 
supporting information). In the limit of low mode num- 
bers p, we find the expected Rousean scaling (X^) oc p~^, 



again demonstrating that the chain structure remains rel- 
atively ideal upon the addition of IL. Thus, no significant 
swelling of the chains can be observed, and the structural 
properties are similar as in neat PEO. On a local scale 
however, the PEO chains become rather contracted due 
to the helical coordination sphere of the lithium ions. 



IV. DYNAMICAL PROPERTIES 

The MSDs of the lithium ions in the individual elec- 
trolytes are shown in Figure |H For P20S • x IL, one in- 
deed observes an increase of the lithium mobility upon 
the addition of IL (Figure 4(a) I. In particular, this in- 



crease becomes visible at 1 — 10 ns, whereas the onset to 
diffusive behavior occurs on the same time scale for all 
systems. Nevertheless, the MSD is much smaller than for 
pure PqS • 2.815 IL, and the crossover to diffusion occurs 
much later. This is a first indication that most likely the 
lithium transport mechanism does not change too drasti- 
cally, which is consistent with the observation that nearly 
all cations coordinate to the PEO chains (Figure ^ . In- 
terestingly, for P2o-axS • X IL, no significant difi^erence 
in the lithium MSD can be found (Figure [4 (b)[ ). 

Figure [5] shows the distribution of the life times p{t) 
of complexes between Li^ and distinct EOs, PEO chains 
and TFSI anions. One observes that both the complexes 
between lithium ions and EOs as well as TFSI molecules 
reach a maximum life-span of about 1 ns. It may seem 
surprising that a significant fraction of the TFSI anions 
is coordinated to Li^ for a few hundred picoseconds, es- 
pecially when keeping in mind that TFSI is a weakly 
coordinating anion and also Figure [3] suggests that this 
coordination is rather rare. One additionally finds that 
the MSDs of Li+ and TFSI have a similar magnitude 
on this time scale, indicating correlated motion of both 
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Figure 4: MSDs of the lithium ions in the (a) P20S ■ x IL and (b) P20-Q1S ■ x IL electrolytes. The curve for the pure IL-salt 
solution PoS • 2.815 IL is also shown for comparison. 




Figure 5: Distributions p(r) of life times r for complexes between Li"'' and EOs, TFSI anions, or PEO chains for (a) P20S ■ x IL 
and (b) P2o-Qa;S ■ x IL. 



species. However, a previous analysis on a PEO/LiI elec- 
trolyte^^ showed that in this case the anion-cation coor- 
dinations lasted up to several ten nanoseconds, and also 
the respective MSDs had similar values during this time 
interval. Thus, the overall tendency towards Li^-TFSI 
complex formation can still be considered as weak in case 
of TFSI-based electrolytes. 

For the EO coordinations, small r-values of a few 
picoseconds are most probable, and pe.o{t) decreases 
monotonously for larger life times (Figure [SJ. This trend 
arises from both the diffusion along the chain and the in- 
terchain transfer. However, a more detailed analysis re- 
lated to the first mechanism reveals that many EOs reen- 
ter the first coordination sphere after a few picoseconds. 
Therefore, the PEO segments will be affected by the pres- 
ence of the ions for a longer period. This is also reflected 
by the long-term decay of peo(''') (i- e. r > 10 ps), which 
can be described by a significantly stretched exponen- 
tial (/? = 0.58 with an average relaxation time of 23 ps). 



Due to the breaking an reformation of EO-Li^ bonds, 
the overall decay becomes slower than according to an 
ideal exponential decay. 

The life times of complexes between Li^ and PEO 
chains show a very wide distribution spreading over sev- 
eral orders of magnitude. In some cases, a PEO chain 
coordinates to a lithium ion only briefiy with one or two 
EOs, whereas in other cases the cation remains attached 
to one particular chain for the entire simulation length 
of 200 ns. 

Finally, one notices that all these characteristics hardly 
depend on the IL concentration. In combination with the 
structural observations, this emphasizes that the typi- 
cal coordinations and the associated exchange processes 
roughly remain the same as compared to the binary elec- 
trolyte P20S. As a result, one might speculate that also 
the lithium transport mechanism does not change signif- 
icantly when incorporating ILs into PEO-based polymer 
electrolytes, at least on a qualitative level. 
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V. EXTRACTION OF THE TIME SCALES 

In the following, we will extract the time scales ri, T2 
and T3 as defined by the transport model, and interpret 
the observed differences in the lithium MSD (Figure |3]) 
via these parameters. 



A. Diffusion along the Chain 



As already discussed in several simulation stud- 
iesii^i^i^, the cations move along the PEO back- 
bone, thereby performing a quasi-one-dimensional ran- 
dom walk. Figure |6] shows the mean squared change 
(An^(t)) of the average EO index n. Within this con- 
text, a distinction was made between cations that were 
bound to one specific chain (irrespective of possible other 
coordinations) and cations bound to two chains during 
the whole observation time t. 

Starting from about 100 ps, the dynamics crosses over 
to a regime that is only slightly subdiffusive with an ex- 
ponent of a « 0.8. Qualitatively, this behavior is found 
for all IL concentrations. For sufficiently long chains and 
a significant amount of PEO— Li"*" complexes that ex- 
ist throughout the entire observation time, {An'^{t)) will 
show diffusive behavior on longer time scales. Naturally, 
in the limit t ^ 00, one would expect a crossover to a 
plateau for finite chain lengths. From FigurelHl neither of 
these two effects can be found, indicating that most life 
times of the PEO— Li"'' complexes are too short (see also 
Figure [SJ. Moreover, the lithium ions move on average 
7 — 8 monomers only during 10 ns. Keeping in mind that 
the lithium ions are typically bound to 4 — 5 monomers 
(Figure [3]), these findings imply that the ions have barely 
left their own coordination sphere during the accessible 
time scale. This is also consistent with the fact that the 
same a- value is observed when calculating (An^(t)) de- 
pending on the starting position of the lithium ion, i. e. in 
the center or at the end of the chain (not shown) . There- 
fore, finite size effects of the PEO chains are irrelevant in 
the present case. 

For P2oS- X IL, no significant change of (An^(t)) can be 
observed when varying the IL concentration, indicating 
that the motion along the chain is basically determined 
by the characteristic motion of the PEO backbone, and 
that the surrounding molecules (i. e. PEO chains or IL 
molecules) have virtually no influence on this mechanism. 
Interestingly, the magnitude of (An^(t)) is essentially the 
same for lithium ions bound to one or to two PEO chains 
(not shown). 

Switching to P2o-qxS • x IL, one observes that the 
motion along the chain becomes slower with increas- 
ing IL concentration, especially for P12S • 1.111 IL and 
PgS • 1.667 IL (Figure [6(b)| ), which can be attributed to 
the lower amount of free EOs in the system, as the PEO 
molecules are successively substituted by IL. As a re- 
sult, each PEO chain coordinates to more cations, and 
one might expect that the attached cations repulse each 



other. Interestingly, the value for a remains basically the 
same in this case, and only the magnitude of the func- 
tional form {An?{t)) cx becomes slightly lower. Thus, 
this decrease is obviously related to a reduced mobility 
of the PEO backbone, since one would observe an earlier 
crossover to a plateau for (An^(t)) and thus a lower a- 
value in case of the correlated motion of several cations 
bound to a single chain. 

Although the motion along the chain is still slightly 
subdiffusive for large time scales, approximate diffusion 
coefficients Di{t*) were extracted from Figure [5] in order 
to determine ri. The Di were obtained from (Ari^(t)) 
from the Einstein relation 



{An\t*)) 
2P 



(1) 



where t* represents the time where the Di were deter- 
mined. Ideally, one would choose t* = T3 (determined 
below) in order to estimate the net effect of this mech- 
anism. Unfortunately, the statistics of (An^(t)) are too 
bad on these time scales. Therefore, Di{t3) was deter- 
mined from Figure |6] by extrapolation. The time scale 
Ti was then calculated according to the formalism of the 
transport model^: 



(2) 



Table |lT] summarizes the resulting n-values. For 
the P2o-Q!a:S • X IL electrolytes, one clearly observes an 
increase of ti with increasing IL concentration, since 
more and more EOs become involved in complexes with 
lithium, and the PEO backbone becomes less mobile. 
Contrarily, for the P20S • x IL electrolytes, ri decreases 
only slightly, and the motion along the chain basically 
remains unaltered by the presence of the IL. However, 
one should bear in mind that the ri-values implicitly de- 
pend on T3, and it is not clear if the extrapolation with 
a sa 0.8 underestimates ti, as (An^(t)) might cross over 
to a diffusive type of motion on larger time scales. Never- 
theless, one can conclude that the experimentally^iii and 
numerically (Figure 14]) observed increase of the lithium 
diffusivity for P2oS- x IL must arise from an enhancement 
of at least one of the other mechanisms. 



B. Segmental Motion 

In order to quantify the segmental polymer motion, the 
MSDs (AR2(t)) of the EOs relative to the center of mass 
of the PEO chain have been calculated (Figure [7]) . This 
quantity has been computed for all EOs (i. e. irrespective 
of the presence of an ion), for EOs bound to a lithium 
ion and for the respective attached ions. The criterion to 
consider a cation or EO as bound was that the average 
EO index of the ion did not change more than one, i. e. 
\An(t)\ < 1 for all time frames during t. For the bound 
EOs, no further distinction between additional coordina- 
tions of the lithium ion to e. g. another PEO chain or a 
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system 


n [ns] 




TR [ns 


T2 [ns] 


Cn [A^ ns-i/^] 


C2 [A^ ns-i/2] 


PEO 




1979 


22 




151.5 




P20S 


147 


1662 


45 


167 


89.0 


46.2 








P2oS- 


X IL 






P20S • 0.66 IL 


140 


1570 


37 


89 


92.7 


59.8 


P20S • 3.24 IL 


127 


1571 


24 


68 


115.2 


68.4 






P 


20-Qa:S 


• X IL 






PieS ■ 0.556 IL 


181 


1479 


43 


145 


81.0 


44.1 


P12S • 1.111 IL 


208 


1359 


35 


145 


82.5 


40.5 


PgS • 1.667 IL 


301 


1151 


28 


104 


78.1 


40.5 



Table II: Parameters characterizing the two intramolecular transport mechanisms (see text for further explanation). The values 
for pure PEO are also shown for comparison. 




10^ 
f [ps] 



Figure 7: MSDs of the average EOs (circles), bound EOs (diamonds), lithium ions bound to these EOs ([An (t)| < 1, solid lines 
in the inset). All MSDs have been computed in the center-of-mass frame of the PEO chain. The dashed and the dotted lines 
show the respective Rouse fits. 
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TFSI molecule was made. Thus, these effects are already 
implicitly contained in the curves in Figure [71 

The average EOs (circles) show typical Rouse-like mo- 
tion with the characteristic relaxation time tr . The dy- 
namics of the bound EOs (diamonds) is qualitatively the 
same but protracted. Therefore, it is possible to char- 
acterize the intramolecular dynamics of the bound EOs 
by using a larger, effective Rouse time T2. The lithium 
ions attached to these EOs (inset, solid lines) closely fol- 
low the bound EOs, which gives evidence for their coop- 
erative motion. On short time scales, the MSD of the 
EOs is larger than the lithium MSD due to the addi- 
tional internal degrees of freedom of the PEO backbone, 
but the MSD of the bound cations catches up on longer 
time scales (approximately 1 ns). Due to the collective 
motion, T2 characterizes the dynamics of the attached 
lithium ions as well. 

Figure [7] also shows the Rouse fits 



(AR2(t)) = 



ml 



\l - exp 



(3) 



for the average (dashed lines) and the bound EOs (dot- 
ted lines). It is important to mention that the precise 
value of T2 obtained from these fits also depends on the 
value of (Rq) entering the prefactor of Eq. [3l In order 
to obtain a fit consistent with the plateau value at large 
i, the MSDs of the average EOs were fitted using two 
parameters, i.e. and (Re)eff. In this way, the val- 
ues in Table U] are replaced by an effective mean squared 
end-to-end vector (R^)off, characterizing the maximum 
accessible intramolecular distance for the cations. Once 
(R^)off was determined from the MSDs of the average 
EOs, the MSDs of the bound EOs were fitted using this 
value in combination with a single fit parameter T2 only. 
The fitting parameters (R^)cff , tr and T2 are summarized 
in Table El 

In case of P20S • x IL, (R^)cff decreases only slightly 
when IL is added. Compared to pure PEO, the (R^}off 
are about 16 — 20 % lower in these systems and show good 
agreement with the respective values in Table HI This 
demonstrates that the PEO chains in these electrolytes 
behave relatively ideal, thus additionally validating our 
analysis in terms of the Rouse model. Both tr and T2 
decrease significantly, clearly indicating that the dynam- 
ics of the PEO segments becomes faster with increasing 
IL concentration. Therefore, the IL can be regarded as 
plasticizer in present case. Due to the enhanced polymer 
mobility, the lithium ions also move faster while they are 
coordinated to a specific chain, and the overall lithium 
MSD increases as observed in Figure ID A similar obser- 
vation has been made experimentally for other plasticiz- 
ers like ethylene/propylene carbonate^^i^ or short PEO 
chains embedded in a high-molecular weight matri x^i^° . 
Naturally, this effect has also been observed for ILs in 
the experiments^ that motivated the present work. 

For P2o-Qa;S • X IL, (Ro)off clearly decreases when the 
PEO chains are successively substituted by IL molecules. 



Consequently, more and more lithium ions coordinate to 
a specific PEO chain, making the backbone more rigid 
and thus diminishing the maximum displacement acces- 
sible via the segmental motions. In combination with the 
decline of (R^)cff, the values for tr and T2 also decrease 
when the IL concentration is raised (Table HI]). However, 
it must be pointed out that the tr- and the T2-values of 
different electrolytes can only be quantitatively compared 
in conjunction with the precise value of (R^)eff. There- 
fore, a direct comparison of tr and T2 is only valid for 
similar (R^)off, as in the case of P20S • x IL. This is also 
highlighted by the observations from Figure |7(b)[ where 
the EOs (both types) rather become slower with the IL 
fraction. Therefore, it is convenient to express the effect 
of both (R^)cff and tr or T2 within a single meaningful 
number. To this purpose, the prefactor of the character- 
istic Rouse-like motion for intermediate time scales (i. e. 
TR/iV^ « t « tr), (AR2(t)) = (2(R2)/7r3/2)yi7^, can 
be calculated, making it possible to directly measure the 
magnitude of the MSDs. Table |lT] displays the constants 
Cr = 2(R2)eff/(^3^R)i/2 and C2 = 2{^l),s / {^^T2Y'^ 
for the individual systems, thereby quantifying the mo- 
bilities of the average and the bound EOs, respectively. 
Here, one indeed observes the same trends as from the 
direct comparison of the MSDs: For P20S • x IL, the seg- 
mental mobility expressed by Cr and C2 clearly increases 
with the IL concentration, whereas it slightly decreases in 
case of the P2o-aa;S ■ x IL electrolytes. For these reasons, 
also the lithium dynamics decreases in P2o-aa:S ■ a; IL in 
contrast to P20S ■ x IL. This behavior can again be ex- 
plained by the reduced flexibility of the PEO backbone, 
since the fraction of bound EOs increases when the PEO 
chains are gradually substituted by IL. This effect coun- 
teracts the plasticizing influence of the IL, 

In both investigated types of electrolytes, the center- 
of-mass dynamics of the PEO chains becomes faster with 
increasing IL concentration. This is not too surprising, 
since the motion of the individual chains is less hindered 
by other chains when the PEO concentration is reduced. 

Within the analysis above, the detailed coordination 
of the lithium ions is not specified further. Of course, 
it is questionable if cations bound to two PEO chains 
show the same dynamical behavior as ions bound to one 
chain only, since the former complexes could be regarded 
as transient crosslinks, which would significantly impede 
the polymer dynamics. A more detailed analysis (see 
supporting information) indeed revealed that there is a 
fundamental difference between these two coordination 
types. While for lithium ions coordinated to a single PEO 
chain ions that diffuse along the backbone are signifi- 
cantly faster than those that remain bound to the same 
monomers, no difference between these to ion categories 
could be found for cations attached to two PEO chains. 
This implies that the cations bound to two PEO chains 
experience no effective transport due to the diffusion 
along the chain. Rather, the PEO chain moves reptation- 
like along its own contour past the cation, which results in 
a non-zero (An^(t)) (Figure [5]), but does not contribute 



10 



to the lithium transport. However, this feature can be 
easily implemented by weighting the effect of ti accord- 
ing to the ratio of ions bound to one PEO chain only (see 
discussion in the supporting information). 

C. Interchain Transfer 

When only considering the diffusion along the chain 
and the PEO dynamics as possible transport mecha- 
nisms, the ion remains confined to a finite volume char- 
acterized by the radius of gyration of the PEO chains. 
Moreover, the diffusion of the center of mass of the poly- 
mer chain becomes irrelevant in the long-chain limit. 
Therefore, the transfer of the cation between different 
PEO chains is necessary for the long-range cation trans- 
port and thus any macroscopically measurable conduc- 
tivity. Here, previous MD simulations have shown that 
the dynamics of the lithium ion is independent from its 
past after such a transfer processes^, leading to Marko- 
vian behavior for time scales larger than t^,, and that 
these processes can thus be regarded as renewal events 
in the spirit of the Dynamic Bond Percolation (DBP) 
model by Ratner et al.—. 

1. Transfer Mechamsm 

We start with the investigation of the detailed mecha- 
nism of the cation transfer process. During most trans- 
fers, the ion is coordinated to two PEO chains simulta- 
neously (i. e. the leaving chain and the entering chain, 
see also Table H]) as a transition state. In some events, 
the cation is also temporarily coordinated by TFSI an- 
ions only and migrates to another PEO chain in this way, 
which will be discussed below. 

Figure |8] shows a histogram of the life times r of the 
intermediate in which a Li^ is coordinated by two PEO 
chains (solid lines) . One clearly observes that the result- 
ing distribution is very broad. Whereas most complexes 
exist only briefly (i. e. a few picoseconds), some exist over 
nearly the whole simulation length of 200 ns. For all elec- 
trolytes, no significant difference can be found for pir), 
suggesting a rather universal mechanism. This is also 
reasonable when keeping in mind that the lithium coor- 
dination sphere containing two PEO chains is relatively 
compact, and that the surrounding molecules thus hardly 
affect the transfer process once the complex has formed. 
Of course, it is questionable if also the brief coordinations 
in Figure[S]can result in a successful cation transfer, since 
the separation of the leaving chain and the simultaneous 
formation of a new, stable coordination sphere on the en- 
tering chain will need a certain time. For this purpose. 
Figure |S] also shows p{t) for complexes which resulted in 
a transfer where the ion remained detached from the old 
chain for at least 100 ps. This time scale is motivated 
by a more detailed analysis (not shown), in which we 
found that brief transfers followed by immediate back- 



jumps to the previous polymer chain occur up to 100 ps. 
Of course, these events cannot be regarded a renewal pro- 
cess in the strict sense, since the lithium dynamics will 
not become uncorrelated to its past after such an event. 
Thus, we used this criterion throughout our analysis to 
identify real renewal processes. From Figure|S]one indeed 
observes that the probability of such a successful transfer 
is lower for short t. These curves indicate that it takes 
on average approximately 50 ps to form an intermediate 
capable to provide a successful transfer. Once such an 
intermediate complex has formed, its decay obviously is 
rather spontaneous as indicated by the roughly exponen- 
tial shape of p{t) for t > 50 ps. Again, the curves are 
nearly identical for all electrolytes. 

In order to study the transfer mechanism in more de- 
tail, we monitored the progress of the cation transfer as 
follows: When a second PEO chain enters the coordina- 
tion sphere of a cation that is already coordinated to a 
PEO molecule, the EO coordination number n at the first 
chain naturally decreases. For brief coordinations of the 
second chain, the crown-ether-like structure of the first 
chain will hardly be perturbed. On the other hand, for 
longer times, one would expect a rather symmetric coor- 
dination by both PEO chains, leading to smaller n. In 
case of a successful cation transfer, the EO coordination 
number decreases even further to zero. The minimum 
coordination number Umin that is reached during r can 
therefore be used to monitor the progress and the suc- 
cess of the cation transfer. Figure IH] shows the probability 
p("-min) to find a specific minimum coordination number 
n-min as well as the probability to find any coordination 
number n on the first PEO chain during r. One clearly 
observes that on average 2 — 3 EOs of each chain coor- 
dinate to the Li^ (solid lines). This coordination is inde- 
pendent from the IL concentration as expected from the 
dense packing of the coordination sphere. For p(nmin), 
one notices that the brief coordinations corresponding to 
JT-min = 4 — 5 bccome less likely with increasing IL concen- 
tration for P2o-Q2:S- X IL. As a consequence, the relative 
probability that a complex with n^i^ = 2 is formed in- 
creases. For P20S • X IL, no such effect can be observed, 
and p(nmin) remains basically unaltered by the presence 
of the IL. This indicates that the increased probability of 
'T-min = 2 for P20-aa;S • X IL is uot related to a concentra- 
tion effect, which would also be present for P20S • x IL. 
Here, rather the lower fraction of free EOs, the electro- 
static repulsion of the coordinating lithium ions and the 
resulting rigidity of the PEO backbone impeding the for- 
mation of the helical structure with n = 4 — 5 leads to a 
larger value for p(nmin = 2). In total, a complete cation 
transfer with rimin — has the same relative probabil- 
ity of about 10 % for all IL concentrations and EO : Li 
ratios. 

It is worth noting that, in contrast to earlier findings 
from MD simulations of a PE0/LiBF4 electrolyte^ no 
significant influence of the anion on the PEO-to-PEO 
transfer could be observed. Here, it was found that the 
BF4 anion assisted the transfer by stabilizing the inter- 
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Figure 8: Histogram of the life times r of the intermediate complexes involving two PEO chains which facilitate the cationic 
transfer. Solid: all events, dashed: successful events only (see text), (a) P20S ■ x IL and (b) P2o-a,a:S ■ x IL. 




Figure 9: Probability p(n) to find a specific EO coordination number n for a cation coordinated to two PEO chains (solid) as 
well as the probability p(nmin) to find a minimum EO coordination number rimin on the first PEO chain during the life time of 
the complex (dashed). 



mediate complex. Within the present analysis, no influ- 
ence of the IL and thus a higher TFSI concentration on 
the mechanistic details of the transfer could be observed. 
Most hkely, the minor importance of the TFSI anions for 
the transfer arises from the suppressed tendency to form 
ion pairs or higher-order ionic clusters as compared to 
BF4. 

2. Renewal Time 

In order to determine rs, we counted the number of 
transfer processes N^-^ from the simulations. As above, 
brief transfers and successive backjumps after less than 



100 ps were excluded, while on the other hand trans- 
fer processes via the anions were also counted, since 
they serve as a renewal process in the same way. Sub- 
sequently, the r3-values were determined according to 
T3 = imax^Li+/^tr, where tmax = 200 ns is the simu- 
lation length and A^Li+ = 27 is the number of lithium 
ions in the simulation box. 

The individual values for are summarized in Ta- 
ble [nil One observes for both the V2o-ax^ ■ x IL and the 
P20S • X IL electrolytes that increases with increasing 
IL concentration. Since the PEO molecules become more 
and more diluted, this can at least partly be explained 
as a concentration effect. In order to study this effect in 
more detail. Table Hill also shows the partial volume of the 
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simulation box occupied by IL relative to the total vol- 
ume filled with PEO and IL, V[L,rei = + Vreo), 
since both species are in principle able to serve as a sol- 
vent. When comparing electrolytes with similar ViL,rci, 
e. g. P20S • 0.66 IL with PieS • 0.556 IL or P20S • 3.24 IL 
with PgS • 1.667 IL, it is obvious that T3 is larger for 
the P20-axS ■ X IL electrolytes. This indicates that not 
only the reduced PEO concentration expressed via Vil,ig1 
leads to an increase of T3, but also the reduced number 
of free EOs, which are necessary for a successful transfer. 

Table IIIII also shows the probability piL that a TFSI- 
supported transfer takes place as well as the average 
squared distance (ARLi)iL the lithium ion covers while it 
is connected to TFSI only. For high IL concentrations, it 
becomes more likely for a lithium ion to migrate through 
the IL-rich regions of the electrolyte. Since this scenario 
is less likely for P20S • x IL, one can conclude that the 
transfer via TFSI in P2o-aa;S • a; IL is mainly caused by 
the crowded PEO chains and the lower fraction of free 
EOs. However, from the percentages to find a lithium 
ion solely coordinated by anions at any time, i.e. 0.4 % 
for P12S • 1.111 IL and 2.5 % for PgS • 1.667 IL, it is clear 
that the residence time in the IL-rich regions is rather 
short. Nevertheless, the value of (AR^Jil indicates that 
the PEO-TFSI-PEO transitions may substantially con- 
tribute to the overall lithium MSD if piL is large enough. 

In total, it was observed that the transfer process is 
similar in all investigated electrolytes, and that a uni- 
versal mechanism holds for the PEO-to-PEO transition. 
Therefore, the critical step is the encounter of a sec- 
ond PEO chain, or, more precisely, a free PEO segment. 
Once the complex has formed, it decays according to a 
universal distribution of life times (Figure |H]), and the 
probability that the cation is completely detached from 
the old chain (nmin = in Figure IH]) is identical for all 
systems. Nevertheless, in the limit of high IL concen- 
trations, the a change of the overall lithium transport 
mechanism becomes visible, since the transfer processes 
facilitated by TFSI become more important, especially 
for P2o-Q2:S • X IL. However, as long as the displacement 
the ion covers in the IL-rich region is sufficiently small, 
these transfers are also captured by the formalism of the 
transport model. 



D. Comparison with Experiments 

Finally, we compare the predictions of the lithium 
transport model to the experimental data taken from 
ref. Of course, a direct comparison of Dl! as ex- 
tracted from Figure |4] is impossible, since the motion of 
the short PEO chains {N = 54) in our simulations sig- 
nificantly contributes to the lithium diffusion, whereas 
most experimental studies focus on the long-chain limit 
(e.g. Mw = 4.000.000 g/mol in ref. 0). However, due 
the explicit A^-dependence of the three time scales^, i. e. 
Ti oc N^, T2 oc N"^ and T3 oc 7V°, and using the fact 
that the center-of-mass motion vanishes for N ^ 00, this 



problem can be easily circumvented. 

For the scaling of T2, entanglement effects may become 
relevant, which would manifest itself by a steeper depen- 
dence of the polymer relaxation time for chain lengths 
larger than the entanglement length Nc- However, if 
Ts < Tc, meaning that the lithium ion leaves the chain 
before the latter begins to reptate or the constraints im- 
posed by the crosslinks become important, the dynami- 
cal contribution to the lithium ions is still Rousean^. In 
the case of PEO, experiments revealed that the entan- 
glement regime sets in from about N sa 75 46,47 ^ Based 
on these observations, one can estimate Tc according to 
Te = rR,(iV = 75) = rR,(75/54)2. For P20S, this leads to 
Tc « 87 ns, which is substantially larger than T3. Also 
in case of the highly plasticized P20S • 3.24 IL one finds 
To ~ 46 ns > T3. Therefore, the lithium ion leaves the 
PEO chain before the tube constraints become notice- 
able, and the transport model can also be applied in the 
limit of long chains. 

Due to the renewal events (i. e. the interchain transfer), 
Du can be expressed as 



^Li(iV) 



(.gi2(A^,T3))M3 
6T3 



+ DpEo{N), 



(4) 



where (5i2(t3))m3 corresponds to an averaged, elemen- 
tary step length between two successive renewal events 
(abbreviated as M3), and is characterized by the in- 
tramolecular transport mechanisms (i.e. Ml and M2). 
The center-of-mass motion of the PEO chains with length 
N is characterized by the diffusion coefficient UpEO- As 
shown in earlier worki, 512 can be described by a Rouse- 
like expression 



9l2it) 



ml) 



[l-exp(-g) 



p=i 



(5) 



where 1 /ti2 = 1 /ti + 1 /t2 is a combined relaxation rate 
due to both mechanisms. Assuming exponentially dis- 
tributed residence times at the distinct PEO chains, one 
can write down an analytical expression for the average 
MSD at a given chain: 



(51 



(t-3))m3 = — / dtexp ( -— ) gi: 
T3 Jo \ T3J 



it) 



ml) 



00 ^ 



^ T12 



(6) 



Solving this expression numerically, values for Dli were 
obtained for iV = 54 {0^=^"^) and for iV ^ 00 {D^r^°°), 
which are summarized in Table HVl For D^^^^, we ob- 
serve that the values obtained from Eq. [6] are in reason- 
able agreement with the respective diffusion coefficients 
D^^^ directly determined from Figure 21 showing that 
the predictions of our model are consistent with the nu- 
merical data. For P2o-aa:S- X IL, the deviations of D^^^* 
from D^^^ can be attributed to the large uncertainties 
of UpEO, since the PEO chains are removed in this case. 
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system 


r-i [ns] piL [%] 




V[L,rcl 


P20S 




17.1 ±1.3 2.5 


31.2 ±8.5 


0.00 


P20S • X IL 


P20S 


• 0.66 IL 


18.4 ±1.4 1.0 


34.2 ± 12.6 


0.19 


P20S 


■ 3.24 IL 


24.1 ±1.3 8.5 


101.6 ±25.6 


0.54 


P2o-Qi:S • X IL 


PieS 


• 0.556 IL 25.0 ± 2.0 4.6 


33.6 ± 10.6 


0.20 


P12S 


■ 1.111 IL 


28.4 ±2.2 13.7 


79.3 ± 17.0 


0.40 


PgS- 


1.667 IL 


32.7 ± 2.5 38.2 


225.7 ±51.0 


0.60 



Table III: Average residence times rs (the leaving PEO chain had to remain detached for at least 100 ps as a criterion for 
a successful transfer process). Here, piL is the probability that the lithium ion is transferred by the TFSI anions, whereas 
(ARLi)iL corresponds to the average squared distance the ion covers while it is connected to TFSI only. Also shown is the 
fraction ViL.rci of the molar volume occupied by IL. 



system 


£)MSD 






P20S 




3.133 


2.932 


1.947 




P 


20S • X 


IL 




P20S 


• 0.66 IL 


3.698 


3.649 


2.309 


P20S 


■ 3.24 IL 


4.744 


4.313 


2.392 




P20 


— ax S ■ 


X IL 




PieS 


• 0.556 IL 


3.396 


(2.633) 


1.511 


P12S 


■ 1.111 IL 


3.124 


(2.613) 


1.311 


PsS- 


1.667 IL 


2.803 


(2.808) 


1.166 



Table IV: Lithium diffusion coefficient Dli (in A^ns~^) cal- 
culated by the lithium transport model (in particular equa- 
tions [4] and [6)1 . The values in parentheses have large uncer- 
tainties arising from the bad statistics in the diffusive regime 
PEO center-of-mass diffusion (since the PEO chains are re- 
moved), and do not properly reflect same trends as D^^^ . 



In the following, we will compare our results with the 
experimental data {D'^^) at T = 323 K taken from ref.i2|. 
Despite the temperature gap between simulation and ex- 
periment, we observe identical trends: For a similar di- 
lution series, in which the IL is added to P20S, ^Li*^ ^-'^^'^ 
increases significantly with increasing IL content (i. e. 
D^^P = 0.052 A^ns"! for P20S • 1.0 IL, 0.118 A^ns"! for 
P20S ■ 2.0 IL and 0.126 k^ns'^ for P20S • 4.0 IL). In con- 
trast to this, rather decreases when the PEO chains 
are suhstitutedhy IL, and, consequently, the EO : Li ratio 
decreases [D"^^ = 0.070 A^ns"! for PiqS-O.O IL, D"^^ ^ 
0.051 A^ns-i for PiqS • 1.0 IL, D"^^ = 0.063 A^ns^^ for 
PioS • 2.0 IL and = 0-046 A^ns"! for P5S • 1.0 IL). 

As discussed above, these observations can be under- 
stood as follows: In case of P20S • x IL, the IL acts as 
a plasticizer, and the enhanced polymer dynamics leads 
to a faster lithium ion dynamics. On the other hand, 
when the EO : Li ratio simultaneously decreases as for 
P20-aa;S • X IL, the remaining PEO chains become more 
and more rigid, which approximately cancels the plasti- 
cizing effect of the IL. 

Thus, two opposing trends related to the segmental 



mobility can be found in ternary polymer electrolyte- 
ionic liquid mixtures: the slowdown of the polymer mo- 
tion with decreasing EO : Li ratio (also encountered 
in the binary systems) and the plasticizing due to the 
IL, leading to an enhancement of the polymer dynam- 
ics. Naturally, for the technologically relevant long-chain 
limit, the enhancement of the center-of-mass motion of 
the PEO chains is irrelevant. For these reasons, one 
should keep in mind that the polymer chains remain mo- 
bile enough when designing new electrolyte materials. 

Of course, the value of Z^Li consists of essentially three 
contributions (Eq. 21): the intrachain distance ((7i2('''3))m3 
the ion covers while connected to the same chain, the 
center-of-mass motion of the PEO molecules character- 
ized by I?pEO E^nd the renewal events measured by T3, 
which facilitate the diffusive motion of the ions in the 
long-time limit for iV 00. Moreover, these trends may 
also oppose each other as in case of P20S • x IL. For these 
reasons, it is desirable to quantify in how far the modified 
intrachain transport and the altered renewal rates con- 
tribute to the overall change of Dli relative to the pure 
polymer electrolyte P20S. 

Table fVl shows the percentage by which the intramolec- 
ular contributions (((7i2(''3))m3), the center-of-mass mo- 
tion of the PEO chains (1?peo) and the renewal rates 
(tj"^) as well as the overall DLi-values change relative 
to the pure polymer electrolyte P20S for = 54 and 
N ^ 00. 

For P2oS- X IL one clearly observes that the intramolec- 
ular contribution increases with the IL concentration, 
mainly as a result of the increased segmental mobility. 
For finite PEO chains, the enhanced DpEO-value addi- 
tionally increases the lithium motion. In contrast to this, 
the contribution of the renewal processes decreases for 
these electrolytes. However, the plasticizing effect over- 
compensates this decrease, and increases for both 
TV = 54 and N ^ 00. 

Remarkably, the relative intramolecular contribution 
also increases for P2o-Qa:S • x IL, although the enhance- 
ment is weaker than for P20S ■ x IL. At first glance, this 
may contradict the findings presented above, where it was 
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N 


= 54 






iV -> 00 




system 


M12 [%] 


M3 [%] 


c.o.m. [%] 


tot. [%] 


M12 [%] 


M3 [%] 


tot. [%] 


P20S 


±0.0 


±0.0 


±0.0 


±0.0 


±0.0 


±0.0 


±0.0 


P20S • X IL 


P20S ■ 0.66 IL 


+21.7 


-7.1 


+39.9 


+24.5 


±27.6 


-7.1 


±18.6 


P20S ■ 3.24 IL 


+53.9 


-29.0 


+98.7 


±47.1 


±73.1 


-29.0 


±22.9 






P20-ai:S ■ X 


IL 








PieS • 0.556 IL 


+9.3 


-31.6 


+10.2 


-10.2 


±13.5 


-31.6 


-22.4 


Pi2S- 1.111 IL 


+6.4 


-39.8 


+23.1 


-10.9 


±11.8 


-39.8 


-32.7 


PsS • 1.667 IL 


+4.7 


-47.7 


+51.6 


-4.2 


±14.5 


-47.7 


-40.1 



Table V: Changes of the intramolecular contribution {gi2{rz))mz, the renewal times rs, the center-of-mass motion Dpeo and 
the resulting DLi-value relative to P20S for A*' = 54 and N ^ 00. 



observed that both the diffusion along the PEO chain and 
the segmental motion become slower vifhen IL is added 
to the system. However, this mismatch can easily be re- 
solved by the following considerations: First, the fraction 
of lithium ions coordinating to one PEO chain only in- 
creases when going from P20S to PgS • 1.667 IL (TablelH. 
As discussed above (see also supporting information), the 
diffusion along the PEO backbone becomes more efficient 
in this way, since lithium ions coordinated to two polymer 
chains experience no effective transport via this mecha- 
nism. The second, more important effect leading to the 
observations from Table IVl is the fact that (gi2(T3))M3 
is evaluated for different (Eq. IB]). Thus, with increas- 
ing T3 also the distance (gi2('''3))M3 increases, although 
this dependence is weaker than linear, and the overall 
Du diminishes in this case. Consistently, calculating the 
individual values for ((7i2(''3))m3 with a constant value 
of T3 = 17.1 ns as observed for P20S, the intramolecular 
contribution indeed decreases for P2o-q2;S • x IL. The 
only contribution which displays a clear increase is i^pEO- 
Thus, for N = 54, Du diminishes only slightly, whereas 
a significant decrease can be found in case oi N ^ 00. 

So far, we focused only on the high-temperature limit 
which we can address in our simulations. Although sim- 
ulations at low temperature would be too costly, we are 
nevertheless able to make some statements about this 
regime. 

For the scenario of the addition of IL, it was reported 
that the glass-transition temperature Tg significantly de- 
creases with increasing IL concentration- , which is an in- 
dication that the enhanced polymer dynamics also plays 
a crucial role at low temperatures. Thus, one would ex- 
pect that the plasticizing effect and thus the increase of 
Du becomes even more pronounced. Indeed, at 323 K, 
Du increases by a factor of 2.4 when going from P20S to 
P20S • 4.0 IL2,, whereas at 423 K, the relative increase of 
D^--^°° between P20S and P20S • 3.24 IL is only half as 
large (Table E]). 

In contrast to this, the decrease of Tg is significantly 
weaker between P20S and PioS • 1.0 IL (which slightly 
resembles the case of IL substitution) , and D^^^ approx- 
imately remains constant at 323 K^. In contrast to this. 



one observes from Table IVl that _D^^°° significantly de- 
creases when going from P20S to PgS -1.667 IL. However, 
one has to keep in mind that these values do not contain 
the contribution of the migration through the IL-rich 
regions, which are relevant for the case of substitution 
(Table |IIT|. Thus, the motion of the ions is partially de- 
coupled from the polymer chains, and one cannot solely 
argue on the basis of the Tg-values. This issue might be 
addressed in future work. 



VI. CONCLUSION 

In this work, we applied a previously developed cation 
transport model^ to MD simulation data of two classes 
of ternary mixtures of PEO/LiTFSI and an ionic liquid. 
Motivated by experimental observation o^'^^ , electrolytes 
of the type P2oS- x IL were studied, i. e. in which an ionic 
liquid has been added to the binary polymer electrolyte 
P20S. Furthermore, we focused on the theoretically in- 
teresting P2o-aa:S • X IL elcctrolytcs, in which the PEO 
chains are substituted by the ionic liquid, in order to ex- 
plore the crossover from the binary polymer electrolytes 
to lithium-salt-containing ionic liquids. 

From a structural point of view, the PEO chains them- 
selves were only locally distorted by the crown ether-like 
coordination structure, and no imprint of swelling or col- 
lapsed chains due to the presence of the ionic liquid could 
be found. 

For all electrolytes, it turned out that the lithium dy- 
namics mainly took place at the PEO chains. A transport 
mechanisms in which the lithium ions are decoupled from 
the PEO dynamics due to the progressive coordination of 
TFSI anion becomes relevant only for P12S • 1.111 IL and 
PgS • 1.667 IL, in which the EO : Li ratio was relatively 
low, resulting in a lack of free ether oxygens capable to 
bind the lithium ions. 

The application of the transport model revealed that 
the increased lithium diffusivity in P20S • x IL can be at- 
tributed to the enhanced PEO dynamics, and that the 
ionic liquid thus serves as a plasticizer. As a result, 
also the lithium ions coordinating to the PEO chains be- 
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come faster, which is in agreement with experimentally 
observed diffusion coefficients^. 

For P2o-axS • X IL, no significant increase of the 
lithium diffusivity could be observed, since the plasti- 
cizing effect of the ionic liquid was superimposed with a 
slowdown due to the reduced flexibility of the PEO back- 
bone caused by crowding effects. Moreover, the same 
crowding also leads to an additional increase of be- 
yond the simple dilution as also observed for P20S • x IL. 
Thus, the lithium diffusion coefficient in fact rather re- 
duces with the mole fraction of ionic liquid. 

In summary, one faces the situation that the segmen- 
tal mobility of the PEO chains plays a decisive role in 
ternary polymer electrolytes: Whereas the lithium ions 



slow down the PEO chains, the ionic liquid accelerates 
the polymer motion. Thus, for the design of new battery 
materials, one should attempt to render the latter effect 
the dominating one. 
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VII. SUPPORTING INFORMATION 

A. Radial Distribution Functions 

B. Static Rouse Mode Amplitudes 

C. Influence of the Lithium Ion Coordination 
Sphere on the Segmental Dynamics 

Within the analysis from Figure [7l the detailed coor- 
dination of the lithium ions is not specified further. One 
would expect that the rather brief coordinations of the 
TFSI anions do not change the mechanistic details of the 
polymer-supported lithium transport, and that this ef- 
fect is implicitly contained in T2. On the other hand, 
it is questionable if cations bound to two PEO chains 
even qualitatively show the same behavior as ions bound 
to one chain only, since the former complexes could be 
regarded as transient crosslinks. A more detailed analy- 
sis indeed reveals that there is a fundamental difference 
between these two coordination types. This can be seen 
in Figure [T^l which shows the lithium MSD of lithium 
ions bound to one and to two PEO chains for P20S. A 
second distinction was made if the ions diffused along 
the PEO chain (both intramolecular mechanisms) or re- 
mained bound to the same EOs (polymer dynamics only, 
here, the same criterion as above, \An{t)\ < 1, has been 
applied). In case of cations bound to one PEO chain 



only (Figure 12(a) I, ions undergoing both types of in- 
tramolecular transport are clearly faster than those re- 
maining close to the initial position on the PEO chain. 
For lithium ions connected to two PEO chains, no signif- 
icant difference between these two categories can be ob- 



served in the MSD (Figure 12(b) ). This implies that the 



cations bound to two PEO chains experience no effective 
transport due to the diffusion along the chain. Rather, 
the PEO chain moves reptation-like along its own contour 
past the cation, which results in a non-zero (An^(t)), but 
does not contribute to the overall lithium transport. This 
can also be seen from the MSD of the EOs that are ini- 
tially bound to the respective lithium ions, which is also 



shown in Figure [12(b) Here, these EOs are significantly 



faster than the lithium ions during the entire observa- 
tion range. Contrarily, in case of complexes involving a 
single PEO chain only, the MSD of the cations which 
move along the chain becomes equally large as the MSD 
of the initially bound EOs at about 1 ns. Seemingly, 



the lithium ions become even more mobile for larger t in 
accordance with the predictions of the transport model, 
however, one cannot make clear predictions within the 
error bars. From Figure 12(a) one also observes that 
for short times the slipping motion of the polymer chain 
is dominant even for ions bound to one PEO molecule, 
as the MSD of the attached EOs is substantially larger. 
Thus, the libration of the PEO backbone on microscopic 
time scales also leads to slight contributions to (An^(<)), 
which naturally results in a displacement of the EOs in- 
stead of the attached lithium ions. For all other inves- 
tigated electrolytes, the observations from Figure [12] are 
qualitatively the same. 

With respect to the transport model, in particular 
Eq. [5l this effect can be easily captured. Here, the only 
additionally required parameter is the ratio of lithium 
ions coordinated to one or two PEO chains (Table [Tj): 

9i2{t) = ri PEO .912(0 + (1 - riPEo)92{t) (7) 

For the fraction riPEo of cations bound to one PEO 
chain, Eq. [S] remains valid, whereas for ions bound to 
two chains with ti^2PEO ~^ co only T2 is important. 
Within this approximation, the average dynamics of the 
lithium ions (i. e. averaged over ions bound to one or two 
PEO chains while simultaneously allowing the interme- 
diate exchange between these two coordination states) is 
estimated from the structural property ri peo only. Fig- 
ure [l3] shows that this simplistic picture is inde ed valid 
to a good approximation. Here, Figure 13(a) displays 
the MSD of ions bound to one or two PEO molecules 
in P20S as extracted from the simulations. These curves 
have in turn been used to calculate an approximate av- 
erage MSD according to Eq. [7| which is also shown in 
Figure 13(a) In fact, the agreement of 312 (i) with the 



average MSD (i. e. a lithium ion that remained on the 
same chain irrespective of other coordinations) is nearly 
quantitative. Deviations for larger time scales are due 
to bad statistics. In the same spirit. Figure 13(b)[ shows 
the ratio between the approximate 512 (t) and the average 
MSD directly calculated from the simulations for all sys- 
tems. As for P20S, this ratio is close to unity for all other 
electrolytes. These observations highlight that, since the 
crosslinks are temporary and an exchange between both 
coordination types takes place, the long-time behavior 
of the intramolecular dynamics may be estimated by the 
average in Eq. [71 as nicely demonstrated in Figure 1131 
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Figure 10: Upper: radial distribution functions in P20S • x IL for the atom pairs (a) Li"*" — Opeo and (b) Li^ — Otfsi- Lower: 
radial distribution functions in P2o-aa;S ■ x IL for the atom pairs (c) Li+ — Opeo and (d) Li"*" — Otfsi- 



19 



10 
10' 
5 10° 

A 

10-^ 

10-^ 



Ifc^ ' 






1 


: PEO 




: PEO^^LiTFSI 




r PEO^^LiTFSI 0.66 IL 




: PEO^nLiTFSI 3.24 IL 





10 



10 
10' 
S 10° 

A 

10-^ 
10-^ 



PEO 



PEO^^LiTFSI 
PEOjj^LiTFSI 0.556 IL 
PEOj^LiTFSI LUl IL 
PEO„LiTFSI L667IL 



10 



(a) 



(b) 



Figure 11: Rouse mode amplitudes (Xp) in dependence of the mode number p for the two electrolyte classes (a) P20S 



X IL and (b) P2o-cixS 



X IL. The Xp have been calculated from the expression of the discrete Rouse model, i.e. Xp(f) = 
l/2)/iV))R„(f). 




Figure 12: MSDs of the lithium ions bound to (a) one or (b) two PEO chains in P20S. A second distinction was made between 
cations that diffused along the PEO chain (Ml and M2) and ions that remained bound to the same EOs (M2 only). The MSDs 
of the average EOs and the EOs that are initially bound to the respective lithium ions are also shown. All curves have been 
computed in the center-of-mass frame. For the other electrolytes, the scenario is qualitatively the same. 
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Figure 13: (a) MSD of lithium ions bound to one or two PEO chains as well as the MSD of cations which are bound to the 
same PEO molecule during t (irrespective of further coordinations) for P20S. The MSDs of ions bound to one and to two chains 
have been used to compute an approximate average 512 (t) according to Eq. [T] (b) Ratio of 312 (i) relative to the real average 
MSD as extracted from the simulations for the P20S • x IL (top) and the P20-aa;S • x IL electrolytes (bottom). 
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